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Abstract 

Quantum wires with spin-orbit coupling provide a unique opportunity to simultaneously control 
the coupling strength and the screened Coulomb interactions where new exotic phases of matter 
can be explored. Here we report on the observation of an exotic spin-orbit density wave in Pb- 
atomic wires on Si(557) surfaces by mapping out the evolution of the modulated spin-texture at 
various conditions with spin- and angle-resolved photoelectron spectroscopy. The results are inde¬ 
pendently quantified by surface transport measurements. The spin polarization, coherence length, 
spin dephasing rate, and the associated quasiparticle gap decrease simultaneously as the screened 
Coulomb interaction decreases with increasing excess coverage, providing a new mechanism for 
generating and manipulating a spin-orbit entanglement effect via electronic interaction. Despite 
clear evidence of spontaneous spin-rotation symmetry breaking and modulation of spin-momentum 
structure as a function of excess coverage, the average spin-polarization over the Brillouin zone van¬ 
ishes, indicating that time-reversal symmetry is intact as theoretically predicted. 
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The spin degree of freedom of electrons has attracted much attention over the last several 
years and plays a central role in spintronic device structures [Tj, the spin Hall effect M, 
topological insulators [1] , Rashba-type surface or interface states 0E] and others. Typically 
the propagation of electrons is treated as a single particle problem. This simple concept 
fails, however, when electronic correlation effects become important and spin- and orbit- 
effects are entangled like in hard-magnets [7], j= 1/2 Mott-insulators [8j or in new emergent 
topological superconductors mm- Among others, the spin-ordered Au-chain ensemble 
grown on Si(553) [T2HT3] and also magnetic chain structures on insulating or superconducting 
surfaces [EMU are recent illustrations of this entanglement and have attracted a lot of 
research activities due to their full flexibility for manipulation with atomic precision. 
Reminiscent of antiferromagnetism in Cr [18], the problem of spin-orbit coupling (SOC) 
under consideration of electronic interaction has been revisited recently PH (IS| . It was shown 
that inherent spin splitting such as Rashba-type SOC yields Fermi surfaces with nesting 
between opposite helical states causing a so-called spin-orbit density wave (SODW) state 
which cannot be characterized by independent local order parameters. This new emergent 
phase of matter arises without breaking time-reversal symmetry above a critical value for the 
Coulomb potential U and its order parameter, i.e. the energy gap A, is determined by the 
energy scales of both the SOC strength A and the interaction U [1TJ]. Furthermore, a finite 
gap A protects the SODW from spin dephasing against external perturbations like magnetic 
fields, thermal excitation and doping by excess coverage. In case that the interaction energy 
is large compared to the SOC, however, the SODW can be overturned at the expense of a 
spin density wave (SDW). 

Experimentally a SODW can be identified by the difference between charge and spin or¬ 
der, the nesting of spin polarized states and the concomitant opening of a gap, and the 
depolarization of the spin order as a function of Coulomb screening due to the dephasing 
of spin states. The observation of all these signatures in our experimental data will be ad¬ 
dressed subsequently after a brief introduction into our quantum wire system of choice, i.e. 
epitaxially grown Pb-wires on Si(557) [2QTt22j . 

Most remarkably, at the critical concentration of 1.31 monolayers (ML), a one-dimensional 
gap A=20 meV is opened below 78 K so that the Pb layer becomes insulating in the di¬ 
rection perpendicular to the wires while it remains conducting along the chain direction 
[231 23]. As sketched in Fig. [ljy this band gap opening is caused by Pb-induced refacetting 
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of the (557) orientation to (223) oriented facets, thus introducing a new reciprocal lattice 
vector, g= 2nd~ 1 e, which fulfils the (spin-polarized) nesting condition 2kp = g along the 
[112] direction, which is one of the signatures for the formation of a SODW as discussed 
above (d is the interwire distance, kp the Fermi wave vector, e is the unit vector). This 
condition is superimposed on the local 2D character of the densely packed a — x \/3 
phase of Pb on the mini-(111) terraces, which is reflected by closed Fermi circles rather 
than by open Fermi lines (Fig. [lj^). The Fermi nesting in this system leads to the expected 
renormalization of the Fermi surface [2T]. For the sake of simplicity and the fact that the 
gap is comparably small we show here the quasi-Fermi surfaces, i.e. the circularly shaped 
spectroscopic features seen at Ep - 50 meV. This representation takes also the replica bands 
into account which are seen by the MDCs taken necessarily also slightly below Ep. A large 
degree of geometrical order along the wires, coupled with high ID conductance, is evident by 
the 10-fold superperiodicity (s =10xasi) at the critical coverage of 1.31 ML, shown by high 
resolution scanning tunneling microscopy (STM) in Fig. [1^ and also seen by spot splitting 
in low energy electron diffraction (LEED, see below). Strikingly, the system reveals giant 
Rashba-splitting of Afco=0.2 A -1 (Rashba parameter «r= 1.9 eV A) [25] which is by far 
larger than the expected value for chemisorbed Pb monolayers on isotropic semiconducting 
surfaces, e.g. Pb/Ge(lll) [26j. A k 0 corresponds to exactly g/2, i.e. the periodicity of spin 
ordering is twice than that of charge ordering at this critical Pb concentration, so that spins 
on adjacent terraces are coupled in an anti-parallel manner. This spin dependent splitting 
is marked in Figs. [3J),c by different colors. The different periodicities of spin and charge are 
a necessary prerequisite for non-trivial coupling between them. 

Here we show that in this array of atomic wires bound to the (223) oriented facets on Si(557) 
we have the unique possibility to tailor the interwire coupling by adsorption of minute 
amounts of excess coverage 50. The screened Coulomb potential scales with the size of the 
excess coverage as U(SQ) = U(O)-exp(-qSO), where q is the screening parameter [27J. As we 
will show, while the spin-orbit order is preserved (at least) up to a critical additional coverage 
of 0.1 ML, steps are decorated by this excess coverage not randomly, but by formation of 
superstructures whose periodicity depends on excess concentration of Pb [201 EH]- This 
highlights that the mechanisms of spin and charge order are not the same. The coupling 
between charge and spin order is mediated by ID in-gap states and is reduced with increasing 
excess coverage. This leads to spin depolarization as a function of excess coverage. Therefore 
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this mechanism has much more fundamental consequences than the band shift associated 
with standard doping, an observation that has been made also in other metallic chain systems 
on surfaces [29, HDj, however with somewhat different physical scenarios. A collapse of the 
SODW phase is associated by an almost vanishing spin polarization close to S 0 ~0.2 ML, 
where the Pb wires turn into an anisotropic 2D metal. Quite naturally, SODW formation 
and the reduced spin-charge coupling as a function of excess Pb concentration also influence 
other physical properties that are sensitive to both charge and spin, like magnetotransport. 
As expected, the gradual spin depolarization seen in photoemission is closely related to the 
spin dephasing probed by magnetotransport measurements [31]. This relationship can be 
quantitatively described by SODW theory, as we will demonstrate below. 

Results 

The effect of excess coverage. Adsorption of excess coverage on the perfect Pb-chain 
system, as shown in Fig. 0b changes primarily the periodicities of the Si(557) surface per¬ 
pendicular to the wires. We demonstrate by the results collected in Fig. [2] that spin and 
geometric order indeed do not have the same periodicities, although they are connected and 
commensurate to each other by the underlying terrace structure. Moreover, while the charge 
order follows the geometry, the spin order turns out to be robust and unchanged by excess 
coverage. 

The respective charge and spin order can be seen in Fig. & which shows the electronic 
structure perpendicular to the wires. Here the periodicity of the facet structure is seen 
in spectroscopy by Umklapp scattering. From the intense momentum distribution curve 
(MDC) peak structure of the perfect wire ensemble (Fig. 0>) a modulation vector of 0.4 A -1 
can be deduced in agreement with the spot splitting along the [112] direction visible in 
LEED. However, as mentioned before, the MDC signal reveals an additional substructure, 
which is shifted by 0.2 A -1 and has the opposite spin helicity [25]. In Fig. [2^, we marked the 
extrema of the ^-component of the spin polarization by red and blue colored circles. This 
equidistant splitting is a spectroscopic hallmark for a SODW in its ground state |19j . Note 
that we present throughout the paper spin-integrated MDCs measured with high resolution. 
The color-coded spin components shown are deduced from additional MDCs recorded with 
the Mott detector (see Supplementary Note 1 and Supplementary Figure 1). In order to 
reveal good signal/noise ratios in reasonable time the spin distribution has been measured 
50 meV below the Fermi level. However, the excellent agreement with results from DC- and 
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magnetotransport demonstrate that this limitation does not affect our conclusions at all. 
With increasing 50, we see that the intensities in the MDCs are slightly altered due to 
changing photoelectron scattering conditions related to the geometrical superperiodicities 
(see below). However, and more importantly, the positions of the maxima of spin polarization 
remain unchanged in k-space. In contrast to the geometric order, the original spin periodicity 
found at 50=0 ML of 0.2 and 0.4 A -1 , respectively, prevails until we reach the limit of second 
Pb layer growth at an excess coverage of 0.2 ML. At this excess coverage also the original 
SODW phase collapses, as seen by negligible spin polarization of all three components. Since 
the enhanced spin-texture along the Fermi surface (FS) nesting direction only appears below 
0.2 ML excess coverage and below 78 K, it strongly indicates to be an interaction effect, 
rather than a single particle SOC effect. We therefore conclude that the periodicity of spin 
order up to h@=0.2 ML is exclusively determined by the periodicity of the terrace structure 
of (223) facets. This terrace induced spin order is reached with maximum spin polarization 
at the critical concentration of 0=1.31 ML, and is interpreted to represent the ground state 
of the SODW. 

Looking at the geometric order, the diffraction pattern of the perfect phase (50=0 ML, 
Fig. [2^) reveals several characteristic features, e.g. pronounced facet spots along [112] (blue 
arrows), a/3 x a/3 spots (green circle) and intensity streaks at x 2-positions (yellow ar¬ 
rows) originating from the long range ordered (223) facet structure (step density 21 %), 
densely packed Pb-reconstruction on the (lll)-units, and Si-dimerization taking place at 
step sites, respectively. I 11 high resolution LEED studies the a/3 x \/3 spots show split¬ 
ting of 27t/s=0.16 A -1 due to formation of domain walls giving rise to 10-fold periodicity 
(s=38.4 A) along the wires [20,;2B], as confirmed by our STM study (Fig. [1^,, Supplementary 
Note 2, Supplementary Figures 4,5). Similar observations were made for Pb/Si(553) [32], 
Upon deposition of excess coverage, superstructure spots appear within the initial diffraction 
pattern of the 1.31 ML phase whose periodicity is reduced by integer multiples of the terrace 
lattice constant as a function of 50 (see diffraction line profiles perpendicular to the steps 
plotted in Fig. [ 2 ) 3 ). As an example, for 50=0.1 ML the periodicity along the [112] direction 
is doubled (see red arrows) by decoration of every second terrace [2H]. 

This variation of geometric order induced by excess coverage is coupled with the reduc¬ 
tion of the band gap, as measured previously [27]. The results, reproduced in Fig. g 
show that subbands are formed within the band gap (see schematic of Fig. [2jd) , which 
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gradually fill the gap as a function of <50. The gradual decrease of the band gap follows 
A = 20 meV ■ exp(— 50/0.038 ML), which corresponds to q— 26.3 ML -1 . Regarding the 
influence of geometric periodicities on the spin order observed, there is obviously no direct 
impact of the modified chain structures, but the (223) structure is still observed in LEED at 
all excess coverages and at low temperature (cf. Fig. [2^). This is the only geometric feature 
remaining constant upon adsorption of excess coverage. Therefore coupling between spins 
on adjacent terraces can only be mediated via the band gap. Consequently, the reduction 
of the band gap also reduces the effective coupling between spin order on different terraces. 
In other words, decoration of the step edges with various periodicities and excess cover¬ 
ages increases the screened Coulomb interaction, which couples to spin order via spin-orbit 
interaction. This is exactly the SODW scenario. 

Spin dephasing in ARPES. Whether this concept holds can be tested by looking at 
excitations of the ground state represented by the spin order at the critical Pb concentration 
of 1.31 ML. Regardless of the protection of the SODW phase the spin polarization of the 
surface bands should be affected by the excess coverage since the spin dephasing time is 
proportional to the size of the electronic gap [19] • Fig. [3^, shows ^-polarization curves 
obtained along the [112] direction, revealing a pronounced spin polarization along the wires 
([110] direction). Up to <50=0.1 ML the polarization curves reveal an harmonic but damped 
oscillatory behavior. The period of 0.4 as well as the almost symmetric amplitudes, is 
giving rise to a vanishing net-polarization within the Brillouin zone (P=0), indicating the 
presence of time-reversal symmetry. This observation excludes the possibility of a SDW and 
further favors the SODW scenario [19], 

In order to quantify the spin depolarization we have plotted the peak-to-peak values of 
the most intense oscillations for all three spin components (Fig. [3 |d). The ^-component de¬ 
creases exponentially with increasing excess coverage and can be well described by S x (5 O) = 
0.56 • exp(—<50/0.036 ML) + 0.12, i.e. q= 27.8 ML -1 . The (solid) line is a fit according 
to the SODW theory (see below). We want to emphasize that the gradual decrease of the 
band gap A deduced independently from DC-transport measurements reveals an almost 
identical g-value. As shown in Fig. [3jz> the S y - and S z - components are considerably weaker 
(see also Supplementary Figure 2). For an ideal Rashba system these components should 
vanish, while in presence of excess coverage this rule is lifted and might be indicative of 
spin spirals [11]. After the band gap has vanished a fast dephasing of the spin-polarized 
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photoholes within the SODW state mimics almost unpolarized surface bands. As expected, 
the effect of depolarization comes along with a reduction of the coherence length £. Fig. [ 3 J 3 
shows £ of the photoelectrons deduced from the full width of half maximum FWHM=2£ _1 
for two MDC peaks as a function of excess coverage which follow a similar trend as the spin 
polarization. 

Moreover, the spin texture along the wires for the quasi-perfect structure has also been 
measured. As obvious from the MDC shown in Fig. |4}r, the most intense features stem from 
Fermi circles around T-points of adjacent Brillouin zones and their replicas as depicted by 
the corresponding labels in Fig. |4]x In agreement with the analysis of transport presented 
below only the central part of the momentum space is relevant. The 2n/s -modulation 
reflecting the domain wall configuration along the wires (cf. Ref. [23]) is seen mainly by 
further satellite peaks shown in Fig. |4jo. Only their consideration allows a comprehensive 
description of the spin polarizations and corresponding Mott-MDCs (Fig. |4^i). At least for 
the S x - and ^-components significant spin signals are found. Most importantly, the spin 
signal around k ~ 0 A -1 reveals only a single spin component giving rise to a spin polarized 
transport channel as we will show below (^-component shown in Supplementary Figure 3). 
Surface magnetotransport. The dephasing behavior can also be derived from magneto 
conductance, mainly stemming from the electronic states at Ep EH- The scattering times 
for various amounts of excess coverage have been determined from a Hikarni analysis [33] of 
the magnetotransport data and are plotted in Fig. [djl (for further details see Supplementary 
Note 3 and Supplementary Figure 6). In context of the depolarization presented above, the 
spin-orbit scattering time r so is of particular interest. As obvious from Fig. [4j;l, the spin-orbit 
scattering rate deduced for the [110] direction (x, along the wires) depends crucially on the 
amount of excess coverage. The strong suppression of spin-orbit scattering in absence of 
any excess coverage is reminiscent of spin-polarized transport along the wires. Furthermore, 
the spin-orbit scattering rate is reduced by at least 2 orders of magnitude upon deposition 
of excess coverage and reveals a close correlation with the depolarization and decoherence 
behavior seen in angle resolved photoemission spectroscopy (ARPES). We want to point out 
that this is per se not self-evident because ARPES and transport probe electrons at different 
energies. Regardless, it is remarkable that the spin polarization follows the same trend which 
in turn demonstrates the close entanglement of both quantities and is a hallmark of strong 
electronic correlation. As we will show below, both can be consistently explained in terms 
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of the SODW formalism. 

Discussion 

Based on our findings and the fact that within first order the spin lifetime r s is given by the 
spin-orbit scattering time r so in a strongly spin-orbit coupled system, the spin dynamics can 
be described by a kinetic equation [ 55] : 

<5S k , c o < s k > ~S k „ m 

— + S k xf2 k =-+ P 1 

01 T] 

where S k is the spin polarization vector. P is any external spin source, and 7) is the scattering 
time. f2 k is the effective Larrnor frequency defined in our case by f2 k = (D. R k x /kp, 0, Osodw)i 
where Or = 2ankp/h, and Osodw = 2 A /K are the frequencies corresponding to Rashba-type 
SOC and SODW, respectively. In our case, P=0, and the initial condition for the spin is 
jSfc(0) x (the other components of the spin polarization are negligibly small), so we can solve 
the above equation analytically in the limit of Dr < Osodw the Fermi momentum, and get 
S x (t) = S x ( 0) • exp(—t/r so ), where the spin dephasing time is (see for details Supplementary 
Note 4) 

1 + \/l + 0g ODW ?7 2 

T s o = --•" C ° nSt ( 2 ) 

Inserting the exponential dependence of the SODW gap on the excess coverage and keeping 
all other parameters constant, we get a good fit of the spin dephasing time r so to the 
experimental value for a reasonable parameter set of SODW gap A=20 meV, rj =1 xlO 12 s 
and 0r=2.6 xlO 12 Hz which is small compared to OgoDW up to <50=0.1 ML. Using the value 
of q<r= 1.9 eV A as the Rashba parameter found in our previous study [25], this refers to 
an extremely small Fermi-wave vector component along the wires (k^, x ~10 -4 A -1 ). This in 
turn explains the insensitivity of the propagating electrons along the wires against atomic 
sized defects [23i| . 

Focusing on the direction along the wires, the emergent modulated spin texture in real space 
can be described by dynamic spin spirals which are antiferromagnetically coupled between 
the wires as illustrated in Fig. [4^. The antiferromagnetic coupling along the [112] direction 
is a consequence of the helical nesting found in ARPES [25|. The spin texture along the 
direction of wires represents the dynamic analogue to static one-dimensional skyrmions 
as observed by STM in magnetic structures [Ml E23- Our results are representative for 
systems where SOC and electron correlations are of comparable strength and with a negative 






exchange coupling. While the SODW is robust and insensitive to these couplings the spin 
polarization is strongly affected and only in the pure SODW phase highly spin-polarized 
transport can be expected. In this respect, the continuous improvement in the field of ultra- 
high spin-resolved spectroscopy will allow us to investigate this new quantum phase and 
the associated band gap at soonest in more detail [38, 39]. Our investigations have further 
shown that the electrons within the SODW phase are strongly correlated with electrons 
residing as gap states. With a large screened interaction by reducing the excess coverage, 
the electronic states also become insulating, opening up a possibility for the realization of 
SODW order induced Luttinger liquid in quantum wires. 

Methods 

Sample preparation and measurements. Spin- and angle-resolved photoemission 
(SARPES) measurements have been performed at the COPHEE end station at the SIS 
beamline of the Swiss Light Source [3D] - The data was analysed as described in Ref. [41] • 
Atomic wire structures were grown by evaporation of Pb on Si(557) substrates. The Pb- 
evaporator used in this study has been calibrated before with an accuracy of 0.01 ML by 
analyzing numerous so-called Devil’s staircase phases on Si(lll) (1 ML=7.83 x 10 14 atoms 
per cm 2 ) [32]. The morphology of the clean Si(557) substrate as well as the structure of 
the wire ensemble has been checked by low energy electron diffraction (LEED) (see Sup¬ 
plementary Note 2). Details about the preparation and coverage calibration are reported 
elsewhere m ■ All photoemission experiments were performed with p-polarized light with 
a photon energy hu =24 eV at a base pressure of 1 x 10 -8 Pa. By recording intensities and 
spin-induced scattering asymmetries with the two orthogonal Mott detectors for different 
emission angles, spin-resolved momentum distribution curves (MDCs) close to the Fermi 
energy have been measured. In order to reveal good signal/noise ratios in reasonable time 
the spin distribution has been measured 50 meV below the Fermi level. Nearly all ARPES 
measurements have been performed at T —25 K while the excess coverage was deposited 
during the cooling phase at 80 ± 10 K. Further details about ARPES, STM and the surface 
(magneto) transport are reported in the Supplementary Notes 1-3. 
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Figure captions: 

Figure 1: 

Atomic and electronic structure of quantum wire ensemble, (a) Scanning tunneling 
microscopy (+1 V) revealing the interwire spacing d— 1.58 nm of a (223) facet structure 
(scale bar is 2 nm). The period of 5.8 A=\/3 x cos(30°) x a$i correlates with formation of a 
Pb-^3 x v^3 reconstruction (yellow trapeze shows size of the unit cell). The superimposed 
modulation of s=10xasi is seen as a spot splitting in the LEED patterns in Fig. [2] as well, 
(b) Schematic of the band structure in the direction across the wires ([112] direction). 
The two colors denote the spin orientation in each of the subbands. The distance between 
bands with the same spin helicity is g = 2n/cb=QA A -1 . The Fermi nesting driven energy gap 
A=20 meV has been determined by ARPES and transport measurements [231 [22> 33]. Right: 
the equidistant sequence of both spin bands is nicely reflected by the MDC (spin-integrated) 
taken at the valence band maximum along the [112] direction. The spin signature of the 
subbands has been deduced from spin-resolved measurements [25]. (c) The giant SOC 

ensures a constructive interference of Fermi surfaces (FS) with the same helicity in k-space 
along the [112] direction. Renormalization of the FS gives rise to a ’diot-spot” wavevector 
A k 0 . The direction along the wires is discussed in context with Fig. [4}r-c. 


Figure 2: 

Evolution of wire structure with additional coverage, (a) Sequence of the spin- 

integrated MDCs taken along the [112] direction for various excess coverages <50. The 
symbols denote the spin texture and have been deduced from spin-resolved measurements 
(see Supplementary Figure 1). In order to reveal good signal/noise ratios in reasonable 
time the spin distribution has been measured 50 meV below the Fermi level, (b) LEED 
diffraction profiles along [112] direction of the 1.31 ML Pb/Si(557) phase with various h@. 
Formation of long-range ordered chains are marked by arrows, (c) 2D LEED pattern with 
<5©=0 ML and (50=0.1 ML. The latter shows nicely the doubling of the period along the 
[112] direction (red arrows) [28]. Optical LEED patterns are shown in the Supplementary 
Note 2. (d) Schematic showing the appearance of new in-gap states, (e) The gradual 

decrease of the band gap A with increasing (5© deduced from a previous transport measure¬ 
ment m The error bars deduced from the transport experiments are around 10 %. The 
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(red) regression line follows A = 20 meV • exp{— <50/0.038 ML). 


Figure 3: 

Damping of spin polarization and coherence length, (a) Spin polarization along 
the wires measured along the [112] direction for different amounts of excess coverage. The 
spectra are shifted for better visibility (S y ~, ^-components and Mott-MDCs are shown in 
the Supplementary Figure 2). (b) Peak-to-peak max-values of all three polarization vectors 
as a function of excess coverage. The (solid) line is deduced from theory. The dashed 
lines are guide to the eyes. The error bars are deduced from least mean square fits to the 
spin polarized ARPES data, (c) Coherence lengths £=2/FWHM deduced from the FWHM 
(after deconvoluting the spectrometer function) of the two most intense peaks of the Mott- 
MDCs. The data point for the clean system was recalculated from spectra taken at T = 60 K. 


Figure 4: 

Spin structure along the wires and spin-orbit density wave, (a) S x - and S y - 

components along the [llO] direction (S z is shown in Supplementary Figure 3). (b) MDC 
taken along the direction of the wires without excess coverage. The peak around k cs 0 has 
a fixed spin polarization, (c) Quasi Fermi surface including replica bands from all T-points. 
The labels refer to those shown in (a), (b) and Fig. [jj:. The 10-fold periodicity correlates 
with the replicas along the [110] direction, (d) Spin-orbit scattering times as a function of 
excess coverage deduced from magnetotransport measurements [31]. The dashed line is a 
fit according SODW theory for rgo,i, he. along the wires. For reference also the scattering 
time Tso,y for direction across the wires is shown. Inset: log scale, (e) Real space sketch 
of antiferromagnetically coupled spin spirals (^-direction) propagating along the wires (re¬ 
direction). For the sake of simplicity the 10-fold periodicity is not shown. 
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